In an expanding universe, what doesn't expand? 
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The expansion of the universe is often viewed as a uniform stretching of space that would affect 
compact objects, atoms and stars, as well as the separation of galaxies. One usually hears that 
bound systems do not take part in the general expansion, but a much more subtle question is 
whether bound systems expand partially. In this paper, a very definitive answer is given for a very 
simple system: a classical "atom" bound by electrical attraction. With a mathemical description 
appropriate for undergraduate physics majors, we show that this bound system either completely 
follows the cosmological expansion, or — after initial transients — completely ignores it. This "all or 
nothing" behavior can be understood with techniques of junior-level mechanics. Lastly, the simple 
description is shown to be a justifiable approximation of the relativistically correct formulation of 
the problem. 

I. INTRODUCTION 

It is not hard to explain to students that the galaxies are moving apart like pennies glued to the surface of an 
expanding balloon or raisins in an expanding loaf of raisin bread Q, UOJ EiQl The expanding material represents the 
uniform stretching of space. But if space itself is stretching, does this mean that everything in it is stretching? Are 
galaxies growing larger? Are atoms? The usual answer is that "bound" systems do not take part in the cosmological 
expansion. But if space is stretching how can these systems not be at least slightly affected? And what would it mean 
for a bound system to be "slightly affected"? Would the bound system, for example, expand at a slower rate? Would 
less bound systems expand more nearly with the full cosmological rate? 

It turns out that these questions gets a spectrum of different answers from experts caught unprepared. Part of this 
confusion is the indeterminacy of just what the question means|6[. In this article, we will put aside some subtleties, 
we will focus on a clear simple question, and we will find a clear (and interesting) answer. 

The simple question will take the form of a simple model: a classical "atom," with a negative charge of negligible 
mass (the "electron") going around a much more massive oppositely charged "nucleus." We will put this classical 
atom in a homogeneous universe in which expansion is described by an expansion factor a(t), where t is time. Our 
goal will be to find the extent to which the growth of a(t) causes the atom to grow, i.e. , causes the electron orbit to 
increase in radius. 

In the description of the atom, it will be useful to use two sets of spatial coordinates, both of them spherical polar 
coordinates with the massive nucleus at the origin. The first system consists of the "physical" coordinates r, 9, 4> 
in which r is the proper distance from the nucleus to the electron at a given moment of time. The second set is 
"cosmological" coordinates R, 9, <f>; a point at fixed values of R, 9, <fi is a point fixed in the stretching space of the 
universe and taking part in the cosmological expansion. The two coordinate systems are related by 

r = a{t)R. (I) 

The angular coordinates 9 and <j> are the same in both the physical and the cosmological coordinates since we can 
think of the cosmological expansion as proceeding radially outward from the origin. 

The question whether the atom takes part in the cosmological expansion is then the question: "Does the electron 
follow a trajectory of constant r (no atomic expansion), or of constant R (full cosmological expansion of the atom), 
or does the electron do something "in between" ? 

The nature of the expansion is encoded in the functional form of a(t), and the choice of this function is the choice 
of the kinematics of the expanding universe. The question of what does or does not expand is a kinematical question 
fundamentally unrelated to the physics that constrains the form of a(t). For that reason we shall use choices of a(t) 
that lead to the clearest insights, but we will comment on the relationship of these convenient examples to realistic 
expansion laws. 

The answers given by our model contain both expected and unexpected features. An expected feature is that 
the comparative strengths of the expansion and of the electrical binding determine whether the atom expands. An 
unexpected feature is the "all or nothing" effect of expansion. We will see that a sufficiently loosely bound electron 
will expand with the universe; it will move with constant R. A more tightly bound electron will, after some initial 
disturbance of its orbit, ignore the continuing expansion and maintain constant r. Without some contrivance of special 
cases there is no intermediate behavior. 
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This paper analyzes the expanding atom at two different levels. First, in Sec.^the description uses only Newtonian 
mechanics and basic electrostatics, and should be accessible to physics students in the junior year. Expansion effects 
are introduced in this model through a very plausible phenomenological stretching force. This analysis leads to a 
relatively simple differential equation for the orbital radius r(t). With this simple equation it is possible to get an 
especially clear understanding of the "all or nothing" feature of the atomic expansion, 

Second, in Sec. 1 1 1 1 1 the same classical atom is analyzed using the kinematics of general relativity and Maxwell 
electrodynamics in curved spacetime. The result of this analysis is a differential equation for r(t) that differs slightly 
from that in Sec.[n] We show, however, that the difference is not significant. If the atom is chosen to be sufficiently 
small, non-Newtonian effects can be made negligible. Section llVI summarizes the conclusions of the paper. 



II. NEWTONIAN ANALYSIS 



Our model consists of an unmoving massive nucleus fixed at the origin of a spherical polar coordinate system r, <j). 
The position of an electron of mass m orbiting in the equatorial plane 9 — 7r/2 is described by the functions r(t), 
4>(t). Since only radial forces act on the electron, its angular momentum mr 2 d<f> / dt is conserved, and we define the 
constant of motion 

to be the electron angular momentum per unit mass. In the absense of cosmological expansion effects, the equation 
of motion for r(t) is derived in the usual way and takes the familiar form 

— -- = -- (3) 

dt 2 r 3 r 2 

The constant of electrostatic attraction C, in SI units, is Qq/ (ATreom), where Qq is the magnitude of the product of 
the nuclear and electron charges. 

We need now to consider introducing the effect of expansion. According to Eq. QJ, a point fixed in the cosmological 
expansion, i.e. , a point of constant R, 9, </>, has a radial acceleration 

d 2 r d 2 a/dt 2 

171 = r L ■ ( 4 ) 

al expansion a 

It seems plausible, therefore, that we can treat this term as a radial force per unit mass, and add it to Eq. JSJ to 
arrive at 

d 2 r L 2 C d 2 a/dt 2 

-^- — = -— +r 1 • ( 5 ) 

dt z r a r A a 

From the solution of this equation, and the chosen expansion factor a(t), we can find the radial position R(t) of the 
electron by using Eq. (JIJ. If we combine r(t) or R(t) with <f>(t) from the integration of Eq. (0, we arrive at a complete 
description of the orbit in cither physical space or cosmological space. 
We will first apply this equation to the simple expansion law 

a(t) cx t 2 . (6) 

It is convenient to have a = 1 at an initial moment. We could arrange this by taking a to be precisely t 2 and choosing 
the starting moment to be t = 1. But in doing this we would encounter an awkward difficulty because da/dt would 
not initially be zero. 

Ideally we want to watch a more-or-less circular electron orbit and ask whether its physical radius r or its cosmo- 
logical radius R is constant. To have a more-or-less circular electron orbit we must start the motion with more-or-less 
zero radial velocity. But we cannot have both dR/dt and dr/dt initially zero unless we have da/dt initially zero. To 
avoid favoring circular physical orbits (dr/dt initially zero) or circular cosmological orbits (dR/dt initially zero) we 
need to have da/dt — at the start. We also want to stat the motion with both d 2 R/dt 2 = and d 2 r/dt 2 = 0; this 
requires that d 2 a/dt 2 initially vanish. To arrange this we replace Eq. 10 by 

a(t) = l + t 2 tanh(i) . (7) 

In order to have d 2 R/dt 2 = d 2 r/dt 2 — initially, we start the orbit, rather arbitrarily, at r = R = 1, and we choose 
L 2 = C, that is, we choose centrifugal repulsion initially to balance electrical attraction. 
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FIG. 1: Radius as a function of time for a = 1 + f 2 tanh£ for L = 3.46. The figure on the left shows the behavior of the 
cosmological radius R(t); the figure on the right shows the physical radius r(t). For this model, the atom initially shrinks 
relative to the cosmological stretching, then maintains a constant value R ~ 0.0376. The physical radius r monotonically 
increases. 

In Fig^we show the results of computing Eq. © with the expansion law of Eq. (7J, and with L — 3.46. The figure 
on the left shows that the cosmological radius R is asymptotically constant at a fairly small value. The electron orbit 
initially shrinks relative to the stretching of space then, after t « 10 the electron distance to the nucleus increases at 
the same rate as the cosmological expansion. (In this case it is not meaningful to speak of "orbital" motion, since <fi is 
changing negligibly.) The figure on the right tells the same story from the viewpoint of physical space. The physical 
radius r{t) is monotonically increasing. The atom, in effect, has been ripped apart by the cosmological expansion, 
and the electron remains fixed at a point in cosmological space. 




FIG. 2: Radius as a function of time for a = 1+t 2 tanhi for L — 3.5. The figure on the left shows the behavior of the cosmological 
radius R(t); the figure on the right shows the physical radius r(t). For this model, the atom undergoes oscillatory shrinkage 
relative to the cosmological stretching. The physical radius r oscillates around r ~ 1. Initial oscillations (shown for < t < 20) 
are large, indicating large eccentricity in the physical orbital motion. Final oscillations (shown for 9, 990 < t < 10, 000) are 
very small, showing us that the electron is settling into an asymptotically circular orbit. 

We now repeat the model except for an increase of L to 3.5, along with the corresponding increase in C to keep 
the initial centrifugal repultsion balanced by the initial electrostatic attraction. The results, shown in Fig. [2 are 
qualitatively very different from those for L = 3.46. Now the cosmological radius R decreases sharply, and the 
physical radius r remains constant, aside from minor oscillations indicating eccentricity in the orbital motion. The 
bound physical motion shown for L = 3.5 occurs for all values of L larger than some critical value L C rit around 3.46. As 
L approaches L C rit from above the initial eccentricity of the motion increases, but the late-time motion remains that 
of a circular orbit of unit physical radius. There appears to be no intermediate case. A small increase in the strength 
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of the electrical binding discontinuously changes the qualitative nature of the late-time behavior of the motion. 

This "all or nothing" behavior is not an idiosyncracy of the choice of expansion in Eq. Q . Models have also been 
computed with a = 1 + i 3 , and for other choices. In the case a = 1 + t 3 the requirement is automatically satisfied that 
da/dt — d 2 a/dt 2 = initially. For this choice of ait) there is again a discontinuous change in qualitative late-time 
motion with L cr it ~ 5.185. For smaller L (with C = L 2 ) the atom is ripped apart by the cosmological expansion. For 
large L, the late time trajectory is a circle with unit physical radius. 

We can get some insight into this "all or nothing" behavior by considering the expansion factor 

a(t) = e at . (8) 

We will start the motion with dr/dt — and C = L 2 and we will not worry about modifying the expansion so that 
da/dt = initially. It turns out that in this case we get very clear results without such a modification. 
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FIG. 3: Radius as a function of time for a = exp(O.Olt) and L = 0.033301. The figure on the left shows the behavior of the 
cosmological radius R(t); the figure on the right shows the physical radius r(t). For this model, the atom initially shrinks 
relative to the cosmological stretching, then maintains a constant value R ~ 1.32 x 10 -4 . The physical radius r monotonically 
increases. 
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FIG. 4: Radius as a function of time for a = exp(O.Olt) and L = 0.033302. The figure on the left shows the behavior of the 
cosmological radius R(t); the figure on the right shows the physical radius r(t). For this model, the electron orbit is bound in 
physical space, and undergoes eccentric motion around the nucleus. 

First we will establish that again the qualitative late-time behavior depends discontinuously on L. Figure [3] shows 
the same late-time qualitative behavior as Fig.^ After a sharp initial shrinkage, the cosmological distance R between 
the nucleus and the electron stays constant. The physical separation monotonically increases. Thus, after an initial 
transient, the distance increases at the same rate as the cosmological expansion. For the slightly greater L value 
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0.033302, the results shown in Fig. 0]are analogous to those in Fig. [5] The electron moves in a bounded eccentric 
orbit around the nucleus. 

We have established then that we again have a rather sudden change from late-time motion ignoring expansion to 
late-time motion following expansion. Due to a special feature of the exponential expansion we can get a very simple 
insight into why this is so. For a = e at , Eq. (JBJ takes the form 



dt 2 



I? 



C 



= — - + a r . 



(9) 



Since t does not explicitly appear, the equation can be viewed as that for a particle moving in one dimension under 
the influence of an r-dependent potential. More explicitly, the quantity 



1 



L 2 



C 



(10) 



is constant, so 



V 



2r 2 



C 
r 



a 



(11) 



can be viewed as an effective potential. This effective potential, is plotted in Fig. for three values of L; for each 
we set C = L 2 so that there is a balance of centrifugal repulsion and electrial attraction at r = 1 . In Fig. [S] the 
intersection of the r = 1 dashed vertical line with a curve of V gives the value of E for motion, since it is the value 
of V for dr/dt = 0. 
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FIG. 5: Effective potential for a — e at , with a = 0.01, for several values of L. 

These plots make it clear what is happening. For small C = L 2 , or, equivalently, for large a, the potential energy 
curve lacks a peak at intermediate radius. The radial force (the negative gradient of the effective potential) will 
everywhere be outward, and the electron will move outward where the expansion term a 2 r 2 will dominate. For 
smaller C = L 2 , the particle will be "trapped" in the potential well around r = 1. For a — 0.01, the dividing line 
between the two cases is L — L cr ; t ~ 0.03302 for which the maximum at r > 1 is just equal to the value at r = 1. 

This analysis tells us why expansion is "all or nothing" at late times, at least in the case of exponential expansion. 
The potential energy curve either does or does not trap the particle near r = 1. If it does, then the effect of the 
expansion is slightly to change the shape of the orbit. The atom cannot "somewhat" expand at late time. If the 
potential well does not trap the electron, then the electron will go to large r and its dynamics will be dominated by 
the expansion. 

We can generalize the underlying idea in this explanation to a broader but less precise explanation of the "all or 
nothing" effect if a(i) is not an exponential. The cosmological expansion term r(d 2 a/dt 2 )/a increases at large physical 
distances r from the nucleus. The centrifugal and electrical forces decrease. This implies a sort of instability with 
respect to expansion. If the electron moves outward, the expansion term will become more important and push the 
electron further outward. Heuristically, at least, such an instability accounts for the fact that the atom cannot just 
keep expanding "slightly." 



6 



We end this section with a comment on the connection between our examples for a(t), and "realistic" expansion 
laws. The favored expansion law for the present epoch of our universe is probably a oc i 2 / 3 , the expansion law 
that corresponds, in Einstein's theory, to a spatially flat universe filled with matter that has negligible pressureQ. 
Interestingly, for that expansion law a~ 1 d 2 a/dt 2 is negative, and the cosmological expansion term on the right hand 
side of Eq. (J5J is actually a cosmological compression term. This suggests that the electron will always remain in a 
bound orbit, and at late-time will not take part in the cosmological expansion. This has been explicitly verified with 
numerical solutions. 

A second interesting expansion law is the exponential expansion a oc exp at, the expansion law that was so useful 
for insights about the all-or-nothing nature of the late-time behavior. This expansion law corresponds to a spatially 
flat universe that can be interpreted either as dominated by a cosmological constant in Einstein's theory, or a universe 
filled with "dark energy" HQ. 



III. RELATIVISTIC ANALYSIS 



The analysis in the previous section was based on a phenomenological term in Eq. (|SJ) representing the effect of 
expansion. Here we analyze the problem using relativisitic cosmology and Maxwell-Einstein theory. 
We start with a standard form0 for the spacetime metric of a homogeneous isotropic universe 



ds 2 = -dt 2 + a 2 (t) 



■R 2 (d0 2 + sin 2 



1 - kR 2 



(12) 



Here, as in Sec.|H] a(t) is the expansion factor and R is the cosmological radius, with r = a(t)R the physical radius. 
We use units in which c = 1. 

The constant k can be positive, negative, or zero. Its magnitude is a matter of choice, since rescaling k is equivalent 
to rescaling R and hence to rescaling a. It is conventional, therefore, to choose k to be ±1 or zero. In our examples 
in Sec.m however, we have already fixed a to be unity at the start of our observations of the electron orbit, and we 
have fixed r = R = 1 for the initial radius of the orbit. In this case kR 2 represents, aside from factors of order unity, 
R 2 / R 2 where R c is the radius of curvature of a inconstant spatial hypersurface. If R 2 /R 2 is not negligibly small, it 
means that our atom occupies a significant fraction of the universe. Spatial curvature effects related to this are in 
principle completely distinct from those of expansion of the universe. We want our atom to be very small compared 
to the size of the universe, so we omit the kR 2 term in Eq. 1|12[) and set g^R = a 2 (t). 

The first step in the relativistic analysis is to find the correct description of the electrical attraction. For the 
spherically symmetric electromagnetic field of the nucleus there can only be a component F 0R of the electromagnetic 
tensor F^ v . The Maxwell equations F a f = 0, with a = and with a = R give us 

^(^ 2 ), R = = ^(F^ ( , ))t (13) 



so that the solution must have the form 



R 2 a 3 

The R equation of motion of the electron's 4-velocity U a is 



F° R « i • (W) 



U a U* = -U F 0R . (15) 
m 

For motion in the 9 = 7r/2 plane this becomes, after some manipulations, 

dt [ aU dt) a 2 i?3[/o ~ C aR 2 ( lb > 

Here L = is a constant of the motion, and we have chosen the constant C to be analogous to the same symbol in 
Eq. ©. 

At any point along the motion U , the component of the particle 4-velocity, can be evaluated from 



o / l + L 2 /a 2 R 2 

V l-a 2 (dR/dt) 2 ' k ' 
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For a specific choice of a(t) the relativistic solution for the electron orbit in cosmological coordinates R, 4>, simply 
requires solution of the system defined by Eqs. Hl(jfl and l|17|) . From this solution, the orbit in physical coordinates 
r, <p is found with Eq. Q . 

A considerable simplification is possible, however, if we notice that U° is the Lorentz gamma factor for the electron's 
physical speed through cosmological space. That is, U° is l/vl — v 2 where v is the physical speed of the electrons 
with respect to observers comoving with the cosmological expansion, a fact that is somewhat clearer in the relation 



U° 



dR 
~dt 



-a 2 R 2 



dt 



-1/2 



(18) 



than in Eq. (|17JI . The electron motion may be relativistic in the sense that it is affected by expansion, but may 
be nonrelativistic in the sense that its Lorentz gamma factor is never significantly different from unity. The special 
relativistic effects encoded in U° — 1 are interesting in their own right, but are rather beside the point for the focus 
of this paper, the influence of expansion. 

To be sure that the motion is slow, we need to be sure that the initial velocity is small. This is not at all the case in 
some of our examples in Sec.[n] The examples all start with start with a = R = r = 1, and hence[lO| with the initial 
value U = yl + L 2 . In our examples the values of L C rit were not small compared to unity. This appears to tell us 
that special relativistic effects are important in determining the late-time behavior of a model. This appearance is an 
artifact of the simple forms of a(t) chosen in Sec.[n] The forms of a in fact, did not specify the timescale for t. If we 
make the rescaling 



t —> nt 



then Eq. (J5J becomes 



cPr 

dp 



L 



Z2 



L/k C^C/k 2 



C_ 



d 2 a/dt 2 



(19) 



(20) 



As a specific example, if we choose a = 1 + (i/1000) 3 (rather than 1 + i 3 ), and keep the starting value of the orbit at 
r = R = 1, the value of L C rit is approximately 0.005185 (rather than 5.185). For this case the maximum value of for 
L sa L cr j t is U° « 4 x 10~ 5 , and special relativistic effects have a completely negligible influence. 

Put more generally, the guideline for avoiding special relativistic effects is that the electron must start out at a 
nonrelativistic speed. The condition for having L cr i t correspond to a nonrelativistic speed is that the initial size of 
the atom and the characteristic timescale for expansion, satisfy the relationship (size) / timescale<§; is much less than 
the speed of light (which is one in our units). 



IV. CONCLUSION 



We have presented a simple definitive question about the influence of the expansion of the universe on a very 
particular system: a classical "atom." And we have found a simple definitive answer: In the long-term our atom either 
ignores the expansion (after some initial disturbance), and has an asymptotically constant size, or it is disrupted by 
the influence of the cosmological expansion and it stretches in size at the same rate as the universe. The way in which 
the atom responds to cosmological expansion depends on how tightly bound it is by electrical forces. 

In analyzing this problem we have developed a simple description of expansion, that of Eqs. @ and JSJ, that avoids 
both special relativistic effects and effects of spatial curvature. And we have shown that this simple description can 
be justified starting with an a priori correct calculation. This simple description may be the basis of useful exercises 
in an undergraduate mechanics course in which expansion of the universe is not the main concern, but in which it 
can be a useful motivation. 
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FIGURE CAPTIONS 

FIG1: Radius as a function of time for a = 1 + 1 2 tanhi for L = 3.46. The figure on the left shows the behavior of 
the cosmological radius R(t); the figure on the right shows the physical radius r(t). For this model, the atom initially 
shrinks relative to the cosmological stretching, then maintains a constant value R w 0.0376. The physical radius r 
monotonically increases. 

FIG2: Radius as a function of time for a = 1 + t 2 tanht for L = 3.5. The figure on the left shows the behavior 
of the cosmological radius R(t); the figure on the right shows the physical radius r(t). For this model, the atom 
undergoes oscillatory shrinkage relative to the cosmological stretching. The physical radius r oscillates around r « 1. 
Initial oscillations (shown for < t < 20) are large, indicating large eccentricity in the physical orbital motion. 
Final oscillations (shown for 9, 990 < t < 10, 000) are very small, showing us that the electron is settling into an 
asymptotically circular orbit. 

FIG3: Radius as a function of time for a = exp(O.Olt) and L = 0.033301. The figure on the left shows the behavior 
of the cosmological radius R(t); the figure on the right shows the physical radius r(t). For this model, the atom initially 
shrinks relative to the cosmological stretching, then maintains a constant value R w 1.32 x 10~ 4 . The physical radius 
r monotonically increases. 

FIG4: Radius as a function of time for a = exp(O.Olt) and L = 0.033302. The figure on the left shows the behavior 
of the cosmological radius R(t); the figure on the right shows the physical radius r(t). For this model, the electron 
orbit is bound in physical space, and undergoes eccentric motion around the nucleus. 

FIG5: Effective potential for a = e at , with a = 0.01, for several values of L. 



